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We investigate a model of two Kondo impurities coupled via an Ising interaction. Exploiting 
the mapping to a generalized single-impurity Anderson model, we establish that the model has 
a singlet and a (pseudospin) doublet phase separated by a Kosterlitz-Thouless quantum phase 
transition. Based on a strong-coupling analysis and renormalization group arguments, we show 
that at this transition the conductance G through the system either displays a zero-bias anomaly, 
G ~ |V| _2 (v^ _1 ) ^ or takes a universal value, G = ^ cos 2 (^^^j, depending on the experimental 
setup. Close to the Toulouse point of the individual Kondo impurities, the strong-coupling analysis 
allows to obtain the location of the phase boundary analytically. For general model parameters, we 
determine the phase diagram and investigate the thermodynamics using numerical renormalization 
group calculations. In the singlet phase close to the quantum phase transtion, the entropy is 
quenched in two steps: first the two Ising-coupled spins form a magnetic mini-domain which is, 
in a second step, screened by a Kondoesque collective resonance in an effective solitonic Fermi 
sea. In addition, we present a flow equation analysis which provides a different mapping of the 
two-impurity model to a generalized single-impurity Anderson model in terms of fully renormalized 
couplings, which is applicable for the whole range of model parameters. 



I. INTRODUCTION 

Kondo physics plays a fundamental role for the low- 
temperature behavior of a large variety of physical sys- 
tems like magnetic impurities in metals, heavy fermion 
systems, glasses, quantum dots, etc. Its key feature 
is the quenching of the impurity entropy through non- 
perturbative screening by many-particle excitations in 
the associated quantum bath. For magnetic impurities 
in metals this amounts to the formation of the Kondo 
singlet between the localized spin and electron-hole ex- 
citations in the Fermi seai 

Most of the aspects of single-impurity Kondo physics 
are now well understood after theoretical tools have been 
developed that can deal with its intrinsic strong-coupling 
nature.-' 2 '- However, in many physical systems the inter- 
action of different impurities, i.e., multi- impurity Kondo 
physics, is important. For example in heavy fermion 
systems the RKKY interaction between different im- 
purity spins leads to competition between local Kondo 
physics and long-range magnetic order that determines 
their phase diagram^ More recently, related questions 
about coupled two-level systems have gained much in- 
terest in quantum computation, where decoherence due 
to unwanted couplings among qubits and between qubits 
and environment should be avoided; on the other hand 
the intentional coupling of qubits is the key step to per- 
forming quantum logic operations. 

In the present paper we investigate the case of two spin- 
1 [2 Kondo impurities Si , S2 coupled via a Ising coupling, 

H**°* = K,SfSZ. (1) 
The Kondo coupling of each impurity to its bath is given 



by 

H? = 2J ± (Spi 3 + + 2J z S^l 3 (2) 

where j = 1,2 labels the impurity, and Soj is the bath 
spin operator at the respective impurity site. Further- 
more the two baths are disconnected. 

Two impurities coupled both to baths and among each 
other present the simplest realization of the so-called 
cluster Kondo effect, which has been discussed, e.g., in 
context of disordered Kondo-lattice compounds^ Fur- 
thermore, the dynamics of magnetic droplets or domains, 
formed in disordered itinerant systems near a magnetic 
quantum phase transition^ also leads to models of cou- 
pled impurities like the one considered here. Therefore, 
we will also refer to the two coupled impurities as mag- 
netic "mini-domain" . 

In the context of Kondo impurities, an Ising-like cou- 
pling Q can be thought of as an effective impurity inter- 
action for heavy fermion systems with an easy axis. Also, 
Ising coupling appears naturally in quantum dots that 
are coupled via their mutual capacitance^ here the two- 
level systems are pseudospins representing the number of 
electrons on the dots, and therefore SU(2) symmetry is 
broken from the outset A^iiS Equivalently, one can think 
of two two-level systems with transversal coupling, with 
the experimental realization of coupled flux qubits^i We 
will discuss different formulations and applications of our 
model in the body of the paper. 

Coupled impurities or two-level systems have been 
investigated in a number of papers ji2ii&±2ii£ii£iiLi2ii2i 
where most attention has been focussed on the case of 
SU(2)-symmetric direct exchange coupling between the 
impurity spins, K Si • S2. Here, two different regimes are 
possible as function of the inter-impurity exchange K: 
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FIG. 1: Schematic phase diagram for two Ising-coupled 
Kondo impurities 11121 . The vertical axis denotes the ratio 
between single-impurity Kondo temperature (determined 
by the Kondo couplings J z , J± , and the bath bandwidth D) 
and Ising coupling K z . The horizontal axis (J z ) measures the 
anisotropy of the Kondo coupling; in the universal regime, 
Tjp <C D, isotropic Kondo coupling corresponds to J z <C J z r . 
The model has two phases: At small J z and large K z , the 
ground state is doubly degenerate and the two impurity spins 
are locked into a "frozen mini-domain". In contrast, at large 
J z or small K z the ground state is a singlet with local Fermi- 
liquid characteristics. The quantum phase transition is of 
Kosterlitz-Thouless type. In the universal regime, 
implying J± <C D, Tjp is the only low-energy scale of the 
single-impurity problem, and the phase transition occurs at a 
critical K z proportional to Tjjp with the proportionality fac- 
tor depending on J z , as shown in the figure (solid line). The 
critical K z diverges as J z — » J z r . Below the dashed crossover 
line the two spins form an Ising mini-domain: the low-energy 
fluctuations are associated with the pseudospin degree of free- 
dom, i.e., for antiferromagnetic K z the staggered impurity 
susceptibility is much larger than the uniform one. 



for large antiferromagnetic K the impurities combine to 
a singlet, and the interaction with the conduction band 
is weak, whereas for ferromagnetic K the impurity spins 
add up and are Kondo-screened by conduction electrons 
in the low-temperature limit. Notably, there is no quan- 
tum phase transition as K is varied in the generic situa- 
tion without particle-hole symmetry (whereas one finds 
an unstable non-Fermi liquid fixed point in the particle- 
hole symmetric case) ji&AiiSi 

As has been pointed out by Andrei et al.f* the case 
of Ising coupling is different and particularly interesting, 
because for large \K Z \ the two Ising-coupled spins form 
a magnetic mini-domain which still contains an internal 
degree of freedom as the ground state of ifj| mg i s dou- 
bly degenerate (in contrast to the inter-impurity singlet 
mentioned above). For the case of antiferromagnetic K z 



(which we will assume in the following), the two low- 
energy states of the impurities (forming a pseudospin) 
are | ||) and | J.|). As we will show in this paper, the 
fate of this pseudospin degree of freedom depends on the 
strength and asymmetry of the Kondo coupling J be- 
tween the spins and the bath electrons. 



A. Summary of results 

Here, we will summarize our main results which are 
detailed in the body of the paper, and schematically rep- 
resented in the phase diagram of Fig. ^ A brief summary 
of the methods used to obtain these results is given below 
in Sec. EH 

The model of two Ising-coupled impurities, connected 
to two separate fermionic reservoirs (realized, e.g., by 
attaching two separate leads to two quantum dots, 
Fig. EJ), has two ground-state phases associated to ei- 
ther a screened or an unscreened pseudospin. For small 
Kondo couplings J±, J z and large K z , tunneling between 
the two pseudospin configurations, | |l) and 1 is su- 
pressed at low energies, i.e., the mini-domain is "frozen" 
as T — > 0, and the ground state entropy is So = In 2. In 
contrast, for small K z the two impurities are individually 
Kondo screened, resulting in a Fermi-liquid phase with 
vanishing residual entropy. This implies the existence of 
a quantum phase transition for K z ~ where Tjp 

is the single-impurity Kondo temperature. For isotropic 
Kondo coupling, i.e., small J z , this has been previously 
pointed out by Andrei et alw& 

What is the nature of this phase transition? Does 
it occur by breaking up the Ising-coupled mini-domain, 
or rather by strong fluctuations of the preformed pseu- 
dospin? What are the universal properties of this transi- 
tion? The key observation, which helps to answer these 
questions, is that the system can also be tuned towards 
the quantum phase transition by increasing the Ising 
component, J Zl of the coupling of the spins to the en- 
vironment in a regime where K z ^> (see Fig.^). For 
K z > T$ the mini-domain is stable. However, upon in- 
creasing J z a many-particle effect (related to formation of 
a Mahan excitonSi) enhances the tunneling between the 
two configurations, | \\) and | []), of the mini-domain. If 
J z exceeds a critical value, J z r , this tunneling can quench 

the pseudospin even for infinitesimal (equivalent to 
infinitesimal transverse Kondo coupling Jj_). 

The resulting phase of the "fluctuating mini-domain" 
is actually a Fermi liquid with vanishing residual en- 
tropy. Note that the finite-temperature properties in this 
regime are rather different from that of the Fermi liq- 
uid which is obtained for Tjp ^> K z . For large J z and 
small Jj_, the high-temperature In 4 impurity entropy is 
quenched in two stages: first, at the scale T° rj K z , the 
mini-domain "forms" , quenching half of entropy; second 
the strong fluctuations kill the remaining In 2 entropy 
at a much lower scale T* , this scale T* can be identi- 
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fled with a collective Kondo temperature associated to 
pseudospin screening. (Note that this type of two-stage 
screening is completely different from the one occurring 
for two conventional Kondo screening channels with dif- 
ferent strengths In contrast, for Tjp ^> K z the en- 
tropy of the two spins is quenched simultaneously in a 
single step at the scale Tg'. Nevertheless, the two Fermi- 
liquid regimes are adiabatically connected by a smooth 
crossover, which we will show to be identical to the well- 
known crossover in the single-impurity Anderson model 
from the mixed valence into the Kondo regime. 

The quantum phase transition at J z = J^ r turns out 
to be in the Kosterlitz-Thouless universality class. As- 
suming continuity along the phase boundary (which we 
verified numerically), this is true for the transition at ar- 
bitrary J z < J z r . Furthermore, universality immediately 
implies that the "fluctuating mini-domain" regime with 
its characteristic two-stage quenching of the entropy also 
exists for small J z close to the quantum phase transition 
(see Fig.JTJ). 

Physical observables, like the conductance in a quan- 
tum dot setup, show universal behavior in the vicin- 
ity of the phase boundary. Therefore, we can calculate 
them close to J z r (see Fig. QJ, where the phase transi- 
tion takes place in a regime of large K z being accessible 
to a strong-coupling analysis combined with renormal- 
ization group arguments. Depending on the experimetal 
setup, we find, e.g., a conductance anomaly characterized 
by the exponent — 2(\/2 — 1), or a universal conductance 

e 2 cos 2 /(^ 7r ) at the phase transition. We empha- 

size again that these results are valid close to the quan- 
tum phase transition even for small J z where a strong- 
coupling analysis is not possible. 



B. Methods and outline 

To obtain the physical picture and the results described 
above, we use a combination of six different and partly 
complementary methods. 

In Sec. [H] we (i) map our model of Ising-coupled spins 
to a generalized Anderson model by bosonization and 
refermionization techniques. This mapping is used to ob- 
tain the qualitative structure of the phase diagram and 
analytic results for the phase boundary at large J z . Fur- 
thermore, for the generalized Anderson model it is much 
easier to implement (ii) numerical renormalization group 
(NRG), which is presented in Sec. II VI With the help of 
NRG it is possible to determine numerically the phase 
boundaries in regimes not accessible to analytic meth- 
ods. Furthermore, NRG is essential to establish that the 
phase transition at small J z is continuously connected to 
the one at large J z . 

Making use of this adiabatic continuity is the main 
idea of this paper to obtain analytic results for the quan- 
tum phase transtion. By increasing J z we can tune the 
transition from a regime with K z ~ T { k ] to a regime 



with K z > T [ x\ where we can employ (hi) a strong- 
coupling expansion (Sec. 1111(1 . The strong-coupling re- 
sult is analyzed using (iv) perturbative renormalization 
group, or more precisely power counting, taking into ac- 
count the anomalous dimensions created from an orthog- 
onality catastrophe. Using these methods the phase dia- 
gram for large K z ^> Tjp and the precise position of the 
critical point for K z — > oo can be obtained. In addition, 
we can determine the relevant phase shifts and scaling 
dimensions of leading relevant and irrelevant operators 
(Sec. IVIjl . This allows us to analytically calculate the 
conductance and zero-bias anomalies close to the quan- 
tum phase transition for arbitrary values of /K z , see 
Sec. IVIII To analytically obtain the precise shape of the 
phase diagram for large J z , we develop in App. ^ (v) a 
generalization of the Schrieffer- Wolff transformation to 
take into account the short-range density interaction of 
our generalized Anderson model, leading to associated 
power-law singularities. 

In Sec. IVland App. [B] we re-derive some above results 
independently by using (vi) flow equations. The advan- 
tage of this method is that it has a broad range of ap- 
plicability and gives a more natural description in terms 
of renormalized quantities. The flow equation mapping 
nicely establishes the equivalence of the different Fermi- 
liquid regimes of the model. It also allows us to derive 
the full phase diagram analytically for general values of 
the coupling J z that are not accessible with the strong- 
coupling expansion. 

Transport quantities and corresponding possible exper- 
imental realizations are discussed in Sec. IVIII The most 
promising way to implement Ising-coupled (pseudo)spin 
variables is the use of charge degrees of freedom of two 
quantum dots, following Matveev's proposal^, and em- 
ploy a capacitive coupling which takes the role of K z , see 
Sec. lVIlAI 

Most of our methods are based on the mapping of 
the original two-impurity model to a generalized single- 
impurity Anderson model, except for the strong-coupling 
analysis in Sec. 1111 Al which is applied to the original 
model (and thus directly establishes the existence of a 
phase transition). As we will show, the employed map- 
ping provides a particularly clear picture of the under- 
lying physics, e.g., it establishes the universality class of 
the transition, and allows to make further progress us- 
ing flow equations. Thus, the mapping turns out to be 
extremely helpful for obtaining the complete picture pre- 
sented below. 



II. MODEL: VARIATIONS AND 
TRANSFORMATIONS 

In this section we discuss the various formulations of 
the model under consideration, together with the map- 
ping between them, which is based on the well-known re- 
lation between the spin-boson model and the anisotropic 
Kondo modeli 3 i 22 i 23 
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FIG. 2: Schematic plot of a system represented by the Hamil- 
tonian J3J|. Two spins are coupled via an (anisotropic) ex- 
change interaction J to two leads. The spins interact via an 
Ising coupling K z . For an explicit discussion of possible ex- 
perimental setups see Sec. IVHI 

Throughout this paper, we will consider the so-called 
scaling limit where both K z and the single-impurity 
Kondo temperature are much smaller than the high- 
energy cutoff D of the theory, T^,K Z <C D. Keeping 
J z > and D fixed, this implies J± — > 0. Only in this 
scaling limit the models discussed below can be mapped 
upon each other. In general, the position of the phase 
boundary, i.e., the value of the critical coupling K z r de- 
pends on microscopic details. In the scaling limit, how- 
ever, K^ r /T^p depends only on J z which parametrizes 
the renormalization flow in the single- impurity model, see 
Sec. IIV CI This universality is represented in the phase 
diagram in Fig. ^ 

A. Ising-coupled Kondo impurities 

We consider the model 

H K = H? + H? + fff 2 ins , (3) 

where two spins Si and S2 interact by the Ising interac- 
tion JJJi #12 nS - Each of the spins couples to a separate 
fermionic bath, ■, via an anisotropic Kondo Hamilto- 
nian Hf (j = 1, 2), see Fig.0 

Hf = H [9 a j] + £ J n 57*^(0)^* w (0) (4) 

na{3 

where a, (3 are spin indices and Hol^aj] = 
Ek,*£kc{ aj c ka} with V aj (x) = Ek e ~ lkXc k«r The 
exchange coupling is assumed to be the same for both im- 
purities and has an anisotropic form, J„ = (J±, Jj_, J z ). 

The fermionic baths are assumed to be particle-hole 
symmetric with a bandwidth D; for a rectangular band 
the density of states at the Fermi level is then pp = 1/D. 
Our conclusions are not modified by the presence of a 
particle-hole asymmetry. A comprehensive discussion of 
possible modifications of our Hamiltonian ©, e.g., due 
to tunneling between the fermionic baths, and how they 
influence our results, will be given in Sec. I VII 

A model of form |J3J| may be approximately real- 
ized with real spins in the presence of a strong Ising 



anisotropy. In addition, it occurs naturally as a model 
for capacitively coupled quantum dots; 8 i 9 i 10 i 20 where the 
local operators Sf and Sf describe charge states, i.e., 
pseudospin degrees of freedom, on the two dots. Con- 
crete application of our results to such a situation will be 
discussed in Sec. IV11I 

B. Coupled Qubits in ohmic baths 

An alternative starting point for our model of Ising- 
coupled impurities can be formulated in terms of two 
two- level systems (spin-boson models), H SB = Hf B + 
Hl B + Hlf™ with 

Hf* = H [b kj ] + I a* + X k "UK + b h) > ( 5 ) 

fc>0 

where H [bkj] = J2k>o w fc "kj^kj' an< ^ a transversal cou- 
pling between them, 

K, 

trtrans z z z ( a\ 

Here are the bosonic creation operators for heat 
bath #j. A is the bare tunneling matrix element be- 
tween the two levels. The impurity properties are com- 
pletely parametrized by the spectral function J{lo) = 
Sfe — ^k), which we assume to be of Ohmic form, 

J{uj) = 2au}e- UJ / 2uJ ". 

One realization of this model is the interaction of tun- 
neling centers in glasses through higher-order phonon 
exchanged In the context of quantum computation this 
model arises in studies of decoherence of coupled su- 
perconducting qubits: the transversal coupling K z is 
generated through a superconducting flux transporter, 
and the heat baths describe the environment leading to 
decoherence H/S^ Here, the assumption of two different 
baths for the two qubits is justified, e.g., if the baths 
model electromagnetic noise coming from read-out cir- 
cuits, which are separate for each qubit. 

C. Bosonization 

The equivalence of H K and H SB can be explicitly 
shown in the framework of bosonization. Furthermore, 
we will demonstrate that both Hamiltonians can be 
mapped to a generalized Anderson impurity model. We 
will use this mapping extensively, both to solve the mod- 
els numerically within NRG and to identify the position 
and nature of the quantum phase transition in certain 
limits analytically. 

It is well knownSS that both the Kondo model and 
the spin-boson model are equivalent to a generalized 
resonant-level model to be denned below. Our model, 
H K and H SB , however, consists of two coupled Kondo- 
and spin-boson Hamiltonians. The crucial point is that 
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the assumed coupling H^ 116 and _ff*2 ar "% respectively, is 
transformed trivially by switching between these three 
representations. 

We start from the two Kondo Hamiltonians J3J and 
apply the bosonization identity 



1 



: Fa") G 



(7) 



where a is a short distance cutoff, F a j is an anticommut- 
ing Klein factor {{F^, F a ,j,} — 26jj'6 a <j'), and <p a j is the 
corresponding bosonic field with [(j) rT j(x),d x i(j) rT iji(x')] = 
2itiS(x — x')5jji b aa i . Transforming to bosonic charge and 
spin fields, 4> s / c .j — (4>tj ± <t>lj) > the bosonized ver- 
sion of the Kondo Hamiltonians (J2J reads 



Hf = H [<j> cj ]+H [<j> sj ] 



J± 
2ira 



( e -i^(o) s + J? t.^. +h 



^0*^(0) (8) 



where Hq[4>] — vf J ^\(d x ())) 2 assuming a linear disper- 
sion, efc = vpk. The bosonic charge field, <f» c j, decouples 
and is omitted in the following. 

Applying a general Emery-Kivelson transformation^ 



U 7 = cxp ij^2S X j<f>sj(0) 



(9) 



parametrized by 7, the Hamiltonian Hjj transforms into 



Hf = U^HfU*, 



Hf 



H [(ps. 
2ira 



J, 



jv F ] S?d x <j> sj (0) (10) 



Importantly, the Ising coupling Q is not affected by this 
transformation, ifj| ms = U-fH^^U^. 

For two special values of the transformation parame- 
ter 7 the Emery-Kivelson transformation results in par- 
ticularly interesting forms of the Hamiltonian. First 
consider the case when 7 = \/2. The exponents in 
the spin flip term of 1)1 OJl then vanish and Hf can 
be cast into the form of the spin-boson Hamiltonian 
(JSJ). We can now easily identify the coupling constants, 

A = A fc = - V2v F ) , u k = v F k and b qj 

are the Fourier components of d x (j) S j{x) with <j) S j(x) = 

-E _ ' 



fe>0 



2tt f 
feL ^ 



i6^e ifcx )e ka l 2 . The linear 



dispersion and the form of the coupling \ q result in a 
ohmic form of the spectral function J{lo) — 2awe' u / 2Uc 
with a strength 

a=(J z p-l) 2 , (11) 
where p is the density of states, p = 1/(2ttvf). 



The single-impurity Kondo temperature has in general 
a power law dependence on the "tunneling rate" J±, 



oc J]~ 



K ^ u ± 

with a introduced above. 



(for J±_ < J z ) , 



(12) 



D. Generalized Anderson impurity model 

Applying the Emery-Kivelson transformation with 
7 = s/2— 1 results in exponentials in 1)10(1 having the same 
form as in the bosonization identity Q and can therefore 
be expressed as fermions ^j. The refermionized Hamil- 
tonian can be identified with a generalized resonant-level 
models 



H 



RL 



ff [*,-]+V (4^.(0) +h.c.) (13) 
(0)^.(0) : 



+VF ( djdj 



where Sf 



2 ! 



V" = -U= and VF = y/2J z - (a/2- 

V27ra v 

l)/p. '5 and e? are fermionic operators, where ^ repre- 
sents solitonic spin excitations of the original conduction 
band, and d describes the spin degree of freedom of the 
impurity. The coupling W vanishes for J z p = 1 — 1/V2 
(or a = 1/2 for the spin-boson model), the so-called 
Toulouse point of the Kondo model^ in this case Eq. 
(|13|l reduces to the conventional resonant-level model. 
Furthermore, W < for isotropic small Kondo couplings, 
J Z = J±< D. 

In the new variables the Ising interaction takes the 
form K z (d\d 1 — ^){d\d 2 — \)- If we interpret the bath 
index j = 1, 2 as a pseudospin index a =|, j, we can iden- 
tify the total Hamitonian with a generalized single- 
impurity Anderson model 

H A = i/ [* CT ] + F^(4* CT (0)+h.c.)+X z n dT n rfi 

(J 



with rida = d\d a 



In this representation, the Ising 



interaction translates to a local Coulomb repulsion, and 
W corresponds to an interaction of the localized level 
d(j with the surrounding electrons. In the limit K z = 
the Hamiltonian describes the extensively studied x-ray 
threshold problem. 21 On the other hand, at the Toulouse 
point where W — 0, the standard impurity Anderson 
Hamiltonian is recovered. For large K z the d level is 
mainly singly occupied; its "spin" a corresponds precisely 
to the pseudospin degree of freedom of the original mini- 
domain. Note that the particle-hole symmetry of the 
effective Anderson model corresponds to the symmetry 
under a rotation by 7r around the x-axis in spin space for 
the original model. 

The mapping of the original two-impurity model onto 
the generalized Anderson model l|14|l is one of the central 
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results of our paper, and will be extensively used in the 
numerical study of the phase diagram and the interpre- 
tation of the results. 



E. Parameter mapping via phase shifts 

It is important to note that the precise relation of the 
three models H SB , H K , and H A depends on the cutoff 
structure, i.e., on properties at high energies and short 
distances. All formulas quoted above which relate the 
various coupling constants are actually only valid within 
the cutoff scheme underlying bosonization. However, it 
is generally believed that all three models are equivalent 
independent of the cutoff structure, as long as one consid- 
ers only the universal low-energy properties in a regime 
where K z and are much smaller than any other scale. 

We consider now the (non-universal) mapping of model 
parameters within different cutoff schemes. For small 
values of J± (or A and V) it is possible to calculate the 
precise mapping by investigating the perturbation theory 
in J±, A, and V, respectively using the fact that all three 
models map onto a Coulomb gas^S&Sl - we will not at- 
tempt this here because it is difficult to do it analytically 
for an arbitrary cutoff scheme. However, the mapping of 
J z , W, and a can be obtained directly by matching the 
conduction electron phase shifts in the limit J±, V = 0, 
as phase shifts are measurable low-energy properties. 

In the Kondo model @ , we denote the scattering phase 
shift for antiparallel conduction electron and impurity 
spins by 6j z ; for parallel spins the phase shift is then 
— 6j z . Analogously, the phase shift in the resonant-level 
model <|13[) is 5\y if the d- level is unoccupied and — Syy 
if the d-level is occupied. For a clear distinction, here 
and in the following we denote by Pf a density of states 
of a fermionic band with hnite cutoff, whereas p refers 
to a density of states within the cutoff scheme underly- 
ing bosonization. In the latter scheme, the phase shifts 
defined above are directly proportional to the coupling 
constants, Sj^ — nJ z p/2 and Sw = TrWp/2, where the 
density of states p = 1/(2ttWp). If one uses instead a 
model where the high-energy cutoffs arise from a band- 
structure, one obtains 8j = arctan -. ,' ,z ! 2 }l^ ao ^}^ , 

' Jz l-(-Jz/2)Rcgoo(0) ' 

where goo(^) = J2k u-ek+io+ 1S ^ ne l° ca l Green's func- 
tion of the electrons. In case of particle-hole symmetric 
bands, goo(0) = —iirpp, this relation simplifies to 

Sj r = arctan[7rJ z pi?/2]. (15) 

Similarly, W in l|14|l induces for V = a phase shift 
6w = arctan[7rW / pi?/2]. Matching the various models by 
their phase shifts, the relation Wp — \[2J z p — (V2 — 1) 
derived within bosonization, translates into 

S w = arctan[7rW/> F /2] = 6j,V2- |(\/2 - 1). (16) 
Note that this equation is only valid for small J± and V. 



III. STRONG-COUPLING ANALYSIS 

In this section we analyze the behavior of the system 
for small Tjp/K z . After presenting a general argument 
for the existence of a phase transition, we discuss the 
resulting physics in terms of the generalized Anderson 
model (|14|) . Interestingly, two different strong-coupling 
limits emerge which will be described in Sees. Illl'd and 
IIII Dl As detailed below, both strong-coupling limits dis- 
play a phase transition of the Kosterlitz-Thouless type. 
Furthermore, the limits will be shown to commute, and 
the physical regimes are smoothly connected. 

A. Effective Hamiltonian 

To investigate the phase diagram sketched in Fig. we 
consider first the limit of \K g \ — ► oo. We can restrict the 
considerations to K z — > +oo, as results for K z — > — oo 
are similar because the z-component of the total spin is 
conserved separately for the "1" and "2" subsystems. 

In the limit K z — ► +oo the two impurity spins form an 
antifcrromagnetic mini-domain, with configurations | f |) 
and | ||). No fluctuations can occur for K z — oo (or 
Jj_ = 0), therefore the ground state of the full system is 
a doublet. 

We now set up a perturbation theory in the small pa- 
rameter J±/K z , by deriving an effective Hamiltonian in 
the {| tl), | IT)} subspace of the impurities. The lowest 
process connecting the two states | jj,), | IT) is of order 
0(J\/ K z ). Thus the effective Hamiltonian in the strong 
Ising-coupling limit reads 

Hf s =Hf Sfi +H^ (17) 

where H^fe is given by with the perpendicular 

Kondo coupling set to zero, J± = 0. Note that the size of 
J z can be arbitrary. The mini-domain is flipped by the 
term 

= -^(^ + ^*li*Ti* t T2 * i2 + h-c) (18) 

with — J2k c kai- The zero-temperature stability 
of the frozen mini-domain now depends on whether the 
operator is relevant in the renormalization group 

sense. 

Since is comprised of four electron opera- 

tors, its bare (tree-level) scaling dimension is negative, 
dim [i?^ p ] tree = — 1. This might suggest that the doublet 
ground state with residual entropy In 2 is stable. How- 
ever, in the present problem acquires an anoma- 
lous scaling dimension which modifies this conclusion. 
This can be understood as follows: For large J z a flip of 
the mini-domain suddenly changes the phase shifts of all 
electrons in the leads, thus exciting an infinite number 
of particle-hole pairs. This is the well-known orthogo- 
nality catastrophy, 30 leading to an anomalous long-time 
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response of the electrons. In the presence of a sharp 
Fermi edge this results in a so-called x-ray edge singular- 
ity which is reflected in an anomalous scaling dimension 

of h% p og . 

In the following we will determine this scaling dimen- 
sion using Hopfield's rule of thumb, 31 and identify the 
critical J z where the H^g becomes relevant, resulting in 
"quantum-melting" of the frozen mini-domain. 

To adjust the Fermi sea to a new ground state after the 
domain has flipped once, a certain amount of charge An 
has to flow to infinity. Hopficld noticed that collective 
response of a Fermi sea depends in the long-time limit 
only on this An: the corresponding correlation function 
decays as i~( An ) . In our problem, we have to consider 
four different Fermi surfaces (j = 1,2, a =|, J.) each con- 
tributing independently. 

A domain flip is induced by the operator A = 
SfS^l^^l^^iz- According to Hopfield's rule of 
thumb, the correlation function in the absence of domain 
flips is then given by 

(Ai(t)A(0)) H K o ~t- a , (19) 

where 

a= Yl An %- ( 2 °) 

j=l,2;<r=T4 

The transferred charges An^ are easily obtained from 
the phase shifts using Friedel's sum rule. For exam- 
ple, a spin flip induced by A changes the phase of the 
down electrons in bath 2 from Sj z to — Sj z which corre- 
sponds according to Friedel's sum rule to a charge trans- 
fer of 2Sj S! /n. Furthermore, the annihilation operator 
eliminates one charge and the total charge transfer 
in this channel is given by An^ = 28 j J it — 1. Simi- 
lar arguments give Anp = — An^i = An^i = — Anj^- 

Therefore, the exponent a is given by a — 4 (—^- — 1^ • 

This result can also be verified explicitly by bosonization 
following Schottc and Schotte^ 

From Eq. (|19f) . we can directly read off the anomalous 
scaling dimension of the domain-flip Hamiltonian 118|) 
with respect to the "frozen-domain" fixed point: 

dim[I^] = l-f = l- 2 (^-l)\ (21) 

Here, the first term arises from the engineering dimen- 
sion of i? c g P . For small scattering phase shifts dim[i/^ p ] 
is negative, i.e., the domain flip Hamiltonian is irrelevant 
and the doubly degenerate "frozen-domain" fixed point 
described by is stable. Domain flips become rele- 
vant for Sj z > St with 

*r-5(l-^)- <-) 

Beyond St the fluctuations of the mini-domain grow to- 
wards low energies giving rise to a new phase - in this 



regime the pseudospin describing the mini-domain is still 
well defined, but is ultimately screened at low energies. 
The special value of St is well-known as the Toulouse 
point of the single-impurity Kondo Hamiltonian^ 

We have thus established the existence of a phase tran- 
sition in the limit of J±/K z <C 1 which is accessed by 
varying J z . In the bosonization cutoff scheme, the critical 
value is pJ z r = 1 — l/\/2 corresponding to the Toulouse 
point of the individual Kondo impurities in the original 
model ©• 



B. Relation to the generalized single-impurity 
model 

What are the properties of this "fluctuating mini- 
domain" and what is the nature of the quantum phase 
transition separating the two phases? This question can 
be tackled within the generalized Anderson model (11411 . 

According to Eq. I|16|) the coupling W in l|14f> vanishes 
for Sj z — St, i.e., at the Toulouse point. This does, 
however, not necessarily imply that W can be treated 
as a small parameter in the vicinity of the Toulouse 
point. Rather, it is important to realize that two dif- 
ferent strong-coupling limits emerge, depending on the 
order of limits taken in parameter space: 

(i) If we consider the limit W — > at fixed small J±, 
then the critical K z will diverge according to the analysis 
above. In this case, K z is the largest scale, whereas J± 
and V can be treated as perturbations. 

(ii) Alternatively, we can take J± — > at fixed small 
W < 0. Apparently, the critical K z will vanish (although 
K z /Tk 1), and a strong-coupling treatment has to 
consider that W is a large scale. 

Both limits correspond to physical situations. How- 
ever, the scaling limit, discussed at the beginning of 
Sec. [H] and employed in plotting the phase diagram of 
Fig. ^ is reached by taking J± — > at fixed J z and 
T% /K z , corresponding to the case (ii) above. More 
generally, in a three-dimensional phase diagram which 
involves a third axis labelled J± in addition to the ones 
shown in Fig.^ case (i) applies for any finite J± in a cer- 
tain vicinity of the Toulouse point, before the behavior 
crosses over to case (ii) - this crossover scale vanishes in 
the universal limit J± —* 0. 

We will now analyse both cases separately - interest- 
ingly the two limits will turn out to commute. 



C. W -> at fixed V 

Assuming that K z is the largest scale in the problem, 
which corresponds to a strong local Coulomb repulsion 
on the impurity site of the generalized Anderson model, 
we can directly map it onto an anisotropic Kondo model. 
With =: vp* (0)* CT (0) : and S = \d) a <T a pdp, we can 
rewrite W J2v "do-^*<r = WS^n^ - n*;) + -^(n dT + 
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^d|)(^*T + As charge fluctuations are frozen out 

for large K z , the last term can be omitted, + n^i = 1, 
and we obtain 



=H [* a ] 



E 

naj3 



(23) 



with J„ = (Jf, Jf, Jf) where 



off 



W 2 
K z ' 



j: 



cff 



AV 2 
K z 



+ W. 



(24) 



Spin-up and spin-down in this Hamiltonian correspond to 
the two states of the mini-domain. The phase diagram of 
(I23|l is well known: when jf is increased from negative 
values towards zero (keeping J^ ff ^ fixed) one observes 
a quantum phase transition from a ferromagnetic regime 
with In 2 residual entropy to a Fermi- liquid phase where 
the spin is quenched, see Fig. |3J The So = In 2 phase can 
obviously be identified with our "frozen mini-domain" . 
It is stable for jf < - \jf\ or W < W r with 



W c 



8V 2 
K z 



(25) 



For K z — > oo or V — > 0, the phase transition is located 
at W c = or equivalently <5j. = St as anticipated above. 
Eq. (|25|l is the exact result for the phase boundary for 
K z — > oo, provided that V and W are mapped onto Jj_ 
and J z as described in Sec. [H] We note again that the 
limit considered here does not correspond to the universal 
limit J± — > 0, because this would give K z r — ► at any 
fixed W. 



D. V -> at fixed W 

Anticipating that K z r — > in this limit, we need to 
consider a problem where W is a large local energy scale. 
Interestingly, at K z = and V — the four impurity 
states are degenerate even in the presence of W because 
of overall particle-hole symmetry. This degeneracy is 
lifted by K z which (as above) favors the impurity states 
|t|) and | If) (assuming K z > 0). 

We proceed with an analysis similar to the usual 
Schrieffer- Wolff transformation, now in the presence of 
a large W . As usual, hopping processes of second order 
in V produce an effective pseudospin interaction between 
the impurity and the conduction band. However, as the 
impurity occupation in the intermediate states is differ- 
ent from that of initial and final states, the intermediate 
state physics involves an x-ray edge singularity due to 
the sign change of the strong potential scatterer W. 

Analyzing the matrix elements, we find x-ray edge be- 
havior which is cutoff by K z , i.e., the effective generated 
Kondo couplings are of the form V 2 jK\~^ where 



= -2-S w 

7T 



2 A 
— ow 

7T 



(26) 




Frozen mini-domain 



FIG. 3: Schematic renormalization-group flow of the effective 
Kondo model 1231 A— describing the screening of the mini- 
domain pseudospin in the limit of large K z . Here, J| ff = cor- 
responds to the Toulouse point of the original Kondo model. 
The Fermi-liquid fixed point is characterized by pseudospin 
screening; the frozen-mini-domain phase corresponds to a line 
of fixed points with unscreened pseudospin. The thick line de- 
notes the phase boundary, given by Jf = -|Jl ff |. Variation 
of J z around the Toulouse point in the original model (at 
fixed J± and K z ) corresponds to a parameter variation in the 
effective model as shown by the dashed line; the dot is the 
transition point. The transition is in the Kosterlitz-Thouless 
universality class. 



where Sw — nWp/2 is the phase shift from W. 

The detailed derivation of the generalized Schrieffer- 
Wollf transformation is given in App. El Finally, one 
arrives at an effective Kondo model of the form (|23|l , but 
now with couplings 



jcS _ W 2 f ± (S w ) 



j; 



cff 



4F 2 f z (S w ) 



W (27) 



where A is of the order of the band cutoff, and both f± 
and f z are smooth, dimensionless functions of Sw with 



/±(0) = 1 



lim f z (S w ) 



/»(Q) = i. 

(28) 

As above, the model shows a Kosterlitz-Thouless phase 
transition in this effective Kondo model occurs at the line 
Jf = -\jf\,le. 



W 



W 2 {f±(S w ) + f z (S w )} 



(29) 



where j3 is a function of W according to 1 — /3 = 2(1 — a) 
with a — (J z p — l) 2 as introduced above. To obtain an 
explicit relation between W, V, and K z in the vicinity of 
the Toulouse point, we expand in W . Dropping additive 
logarithmic corrections, we have 

ps W , A* 3 ps 1 . (30) 
With this and we can re-write Eq. lj29"|) as 

1 , l-2c 



W 



(31) 
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which is smoothly connected to the result (|25|l obtained 
in the limit W -> of Sec. ITTTCl 

As announced, the critical K z for fixed W depends 

(1) 1 l-2a 

only on T^- y cx V 1 - a A 1 - a . To fix the prefactor, we 

first have to define for an asymmetric Kondo model. 
This is best done using a physical observable like the 
impurity specific heat coefficient 7 = lim^^g C; mp /T of 
the anisotropic single-impurity Kondo model. We employ 



(32) 



where w — 0.4f28 is the Wilson number. 1 At the 
Toulouse point, a = 0, one easily finds 7 = 1/(3/0_fV 2 ). 
Thus, for W — > we obtain: 



(i) 



A' 



urn 2 K cr 



T (i) 

-1.964-^. 

iT5 r 



(33) 



For particle-hole symmetric bands, this can be re- written 
using Eqs. (|T3|) and into: 

T (i) p F t(;7r 2 ^sin 2 (-^) 



.FT? 



g 2 -^-(J?-Jz) 

0.578 p F {J^ r -J z ) (34) 



valid for J z — > J z r . 

Concluding this analysis, we have shown that also in 
the limit V — > (keeping W fixed) the effective Ander- 
son model (|14f) can be mapped onto an effective Kondo 
model, which describes the screening of the mini-domain 
pseudospin. The phase transition is in the Kosterlitz- 
Thouless universality class. The two strong-coupling lim- 
its are adiabatically connected, as the involved impurity 
states and transitions are similar in both cases; in addi- 
tion the equations for the phase boundaries 1(25(1 and 131(1 
match. 

From the mapping to the anisotropic Kondo model we 
can also identify the "fluctuating mini-domain" regime 
with a Fermi-liquid phase. Here, the pseudospin is 
screened below the Kondo temperature, T* , of the ef- 
fective Kondo model l(23|) . Importantly, the Fermi sea 
of the effective model is formed by solitonic spin ex- 
citations of the original model © - this will strongly 
influence the conductance through the system as dis- 
cussed in Sec. IVIII At the Toulouse point, we can es- 
timate T* ~ min(K z ,D)exp[-\K z \/(8V 2 p F )]. Close to 
the quantum phase transition, reached for K z /V 2 — ► 
at the Toulouse point, T* is exponentially small. Sim- 
iliarly, for finite K z r and K z < K z r , one expects for a 
Kosterlitz-Thouless transition the behavios^ 



-b/y/K^-Kz 



where a is a function of and J z \ this form will actu- 



(35) 



ally be used to fit the numerical data of Sec. IIVI 

At the Toulouse line, one can investigate the crossover 
from the "fluctuating mini-domain" regime to the regime 



with Kondo-screened spins at K z = (dashed crossover 
line in Fig.^). As our model is equivalent to an Ander- 
son impurity model 1 (14(1 . this crossover is equivalent to 
the well-known Anderson model crossover from mixed- 
valence to Kondo behavior. This crossover takes place 
at K z ~ V 2 p F , i.e., when t£ ] /K z is of order 0(1). Si- 
miliarly, one finds for J z — > 00 that this crossover is lo- 



cated at T, 



(i) 



K Zl and furthermore for J, 







T]p (as no other 



-k '~ J -L 
one also expects this crossover at K 
low-energy scale exists in this regime). We therefore con- 
clude that the dashed crossover line in Fig. ^ is always 
located at T^ ] /K z - 0(1). 

The schematic RG flow of the effective Kondo model, 
shown in Fig. [21 illustrates the behavior of the Ising- 
coupled two-impurity model in the strong-coupling 
regime. Three observations are important: (i) Both the 
So = In 2 frozen-domain phase and the So = Fermi- 
liquid phase are stable to small perturbations - this con- 
clusion is in agreement with the analysis of Ref. I20L 
(ii) The So = In 2 phase corresponds to a line of RG 
fixed points, (iii) The two phases are separated by a 
Kosterlitz-Thouless transition, characterized by logarith- 
mic rather than power-law behavior of thermodynamic 
observables. 

So far, the conclusions above mainly apply to the vicin- 
ity of the Toulouse point of the single-impurity model, 
i.e., for strongly anisotropic Kondo coupling. In the next 
section we will present numerical results which strongly 
support that the above picture is valid over the whole 
phase diagram. In particular, we shall show that the 
phase transition is of Kosterlitz-Thouless type even in 
the case of isotropic Kondo couplings. Furthermore, no 
other phase transition (than the one indicated in Fig. ^) 
occurs; this implies that the Fermi liquid formed by the 
fluctuating mini-domain is adiabatically connected to the 
Fermi liquid of two individually Kondo-screened spins, 



realized in the limit K z <C T 



(i) 

K ■ 



IV. NUMERICAL RENORMALIZATION 
GROUP ANALYSIS 



In this section, we turn to a numerical investigation of 
the model of Ising-couplcd impurities using the numerical 
renormalization group (NRG) technique^ In principle, 
an investigation of the original two-impurity model @ 
is possible - however, the required two bands of spinful 
fermions are computationally demanding within NRG, 
and results are significantly less accurate compared to 
the NRG treatment of single-band impurity models. 

Therefore, we have decided to study the generalized 
Anderson model of Eq. 1(14(1 . obtained after bosonization 
and refermionization of the original model. Featuring 
only one band of spinful conduction electrons, it allows 
high-accuracy numerical simulations down to lowest en- 
ergy scales and temperatures. 
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A. Parameters 



In the following, we will show numerical results for 
the generalized Anderson model ((14|). for a rectangular 
particle hole symmetric fcrmionic band of width D = 
2\/2, and different values of the hybrization V, density 
interaction W , and on-site repulsion K z - note that K z 
is identical to the original Ising interaction between the 
two impurities, whereas V and W are related to Jj_ 
and J z as detailed in Sec. [H] In particular, W mea- 
sures the deviation from the Toulouse point. According 
to Eq. JUJ, valid for small V, J z — corresponds to 
Wp F = -(2/7r)tan[7r(\/2 - l)/2] w -0.48. Similarly, 
large negative values of W, i.e., W — * — oo, correspond 
to pj z = 1 - V2 = -V2 ( / oJ z )tou1ou SO with pj z here de- 
fined in the bosonization cutoff scheme. 

The mapping between V and Jj_ (which are simply pro- 
portional) can be achieved via the Kosterlitz-Thouless 
transition line for a single impurity: This line, where 
Tjp vanishes, is given by J z = — \J±\. Using NRG, we 
have numerically determined a few points on this line, 
characterized by parameters V and W in the model for- 
mulation (|14fl with K z = - note that this involves an 
extrapolation of Tk(V) to Tk = 0. With the mapping 
between W and pJ z established above, we find for the 
employed parameters the correspondence pj± 0.38V, 
valid for small V. 

Within NRG, the bath density of states is discretized 
on a logarithmic grid, with a discretization parameter A, 
i.e., the energy axis is divided into intervals at points 
D, Dhr 1 , DA~ 2 , . . .. The discretized model is trans- 
formed into a semi-infinite chain form, and then diagonal- 
ized iteratively^ After each diagonalization step the low- 
est N s eigenstates of the Hamiltonian are kept. Clearly, 
the results are "exact" in the limit N s — ► oo, A — ► 1. 
In practice, numerical results are monitored at fixed A, 
where convergence upon increasing the value of N s can 
be readily achieved. Then, these converged numbers are 
extrapolated as function of A to A = 1 to obtain an esti- 
mate for the "exact" result. 

Most of the NRG calculations here have employed a 
discretization parameter of A = 2, keeping N s = 650 states 
per NRG step. For selected parameter values we have 
performed calculations with different A and N s . For each 
A convergence w.r.t. N s can be easily achieved, however, 
for some quantities the A dependence turns out to be 
rather strong. Therefore, the A — > 1 extrapolation has 
to be performed carefully, as will be detailed at the end 
of this Section. 

The results shown below primarily correspond to the 
universal regime of V -C W, D and K z <C W, D; we have 
also performed some calculations in the regime K z 3> W 
(not shown) with results consistent with the analysis of 

sec. nun 



B. Results for RG flow and entropy 

In Fig. 0] we show NRG flow diagrams displaying the 
energies of a few low-lying many-body eigenstates as 
function of the number of NRG steps. The data in Fig.0^, 
clearly show that for small values of K z the same fixed 
point is reached for any V and W - this fixed point 
can be identified with the Fermi-liquid phase, which is 
in particular also reached for K z = 0. Therefore, the 
Fermi-liquid regime of two separately Kondo-screened 
impurities is adiabatically connected to the "fluctuating 
mini-domain" regime which can be characterized by pseu- 
dospin screening below the collective Kondo temperature 
T*. In Fig. flow diagrams for larger values of K z are 
shown - here the fixed points reached at low energies are 
very similar for different parameter sets, but not identi- 
cal - this is consistent with the notion of a line of fixed 
points with In 2 residual entropy. 

It is important to emphasize that no additional fixed 
point is observed for K z « K C J , which could possibly 
correspond to a critical fixed point. This clearly shows 
that the quantum phase transition in our problem is not 
associated with standard critical behavior, but indicates 
that it is of the Kosterlitz-Thouless universality class. 

To characterize the fixed points, we have evaluated the 
impurity entropy S(T) using NRG. In Fig. [5] we show 
results for different values of W and several K z . The 
discussed two-stage quenching of the entropy, occuring 
for Tjp < K z < K z r , can be nicely seen in all panels 
- note that panel c) shows data for W > 0, i.e., on the 
right-hand side of the Toulouse point. 



C. Phase boundary and phase diagram 

From the NRG results for fixed values of W , V and dif- 
ferent K z < K z r it is possible to extract a characteristic 
crossover temperature T* below which the pseudospin is 
screened, see above. For numerical simplicity we defined 
T* through S(T*) = 0.4. The dependence of T* on K z 
allows to determine K z r where T* vanishes. We fitted 
the data with T*(K Z ) = aexp[b/ \/Kf - K z ] © with 
fit parameters a, b, K z r ', which is the form expected near 
a Kosterlitz-Thouless transition^. The fit works excel- 
lently for all anisotropies, as shown in Fig. - this is 
again strong support for the Kosterlitz-Thouless nature 
of the transition. 

The single-impurity Kondo temperature, T^p for given 
W and V is determined from Eq. \\'.V2\\ , where the specific 
heat coefficient 7 is extracted from the NRG data for 
S(T). 

Having determined both tP and K z r , we are in 
the position to plot the phase diagram in the univer- 
sal fashion indicated in Fig. ^ i.e., employing the limit 
tP <C D. The result is shown in the main panel of 
Fig.0 

It is possible to make the meaning of universality more 
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a)5 = 



2 



5 < 



a) Wp F = -0.44 





FIG. 4: NRG flow diagram for the generalized single-impurity 
Anderson model 1141 , for parameter values belonging to a) the 
Fermi-liquid phase with So = 0, b) the frozen-domain phase 
with So = In 2. Solid: Wp F = -0.44, V = 0.075 (Kf = 
7.6xl0~ 10 ), Dash-dot: Wp F = -0.10, V = 1.5xl0~ 3 (Kf = 
3.4 x 10" 6 ), Dashed: Wp F = -0.034, V = 1.5 x 10~ 5 (Kf = 
4.8 x 10 -9 ). In a) and b), K z has been chosen slightly below 
and above the critical value, respectively. For all parameters, 
the system is in a S = In 4 regime at high temperatures (small 
N), in a) it flows to the 5 = state by passing through 
a regime with S = In 2. In a), the additional dotted curves 
show the flow for W p F = -0.44, V = 0.075, and K z = 0. The 
Wp F values span a large range of anisotropics; nevertheless, 
the S = fixed point is unique, and the finite-temperature 
crossover is universal for the curves close to Kf . Panel b) 
nicely shows that S — In 2 actually corresponds to a line of 
fixed points. NRG parameters are A — 2 and N a =650. 



precise: so far we have distinguished the parameter sets 
by their value of J z (or W), leading to different values 



Of T; 



(i) 



and 



■ K '/Kf . Within the RG treatment of Yuval 
AndersonSS for the single-impurity Kondo model, it is 
easily seen that we can expect identical low-energy be- 
havior for two single- impurity models if the initial param- 
eters place the two models on the same RG trajectory. 
(The RG trajectories are identical to the ones shown in 
Fig- El) Therefore, the correct parameter for the horizon- 
tal axis of our phase diagram is a parameter labelling the 
RG trajectories, i.e., a proper RG invariant. Note that 
for J z > and J± — > 0, J z can be used as such a label - 
this is what we have done so far. A proper RG invariant 
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FIG. 5: Temperature evolution of the impurity entropy cal- 
culated by NRG for the generalized single-impurity Anderson 
model 1141 for different anisotropics of the Kondo coupling. 
In the "frozen mini-domain" phase the residual entropy is In 2 
while it vanishes for K z < Kf. For > K z (solid curves ) , 
the high-temperature In 4 entropy is quenched in a single step, 
whereas two-stage screening occurs for < K z < Kf . a) 
Wp F = -0.44, V = 0.15 (Kf = 1.5 xlO" 5 ), close to isotropic 
Kondo coupling. K z is: Solid 0, long-dash 10~ 5 , long-dash- 
dot 1.3 x 10" 5 , short-dash 1.5 x 10" 5 , short-dash-dot 10 -4 . 
b) Wp F = -0.034, V = 1.5 X 10~ 5 (Kf = 4.8 x 10~ 9 ), i.e., 
close to the Toulouse point of the individual Kondo impuri- 
ties. The K z values are: Solid 0, long-dash 10~ 9 , long-dash- 
dot 1.5 x 10" 9 , short-dash 3 x 10~ 9 , short-dash-dot 10~ 7 . c) 
Wp F = 0.44, V = 1.5 X 10~ 7 , i.e., on the right-hand side of 
the phase diagram Fig. Qwhere no phase transition occurs as 
function of K z . K z is: Solid 0, long-dash 10 -8 , long-dash-dot 
10" 7 , short-dash 10" 6 , short-dash-dot 10" 5 . 
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FIG. 6: Numerically determined values of T*(K Z ) for WpF = 
-0.034, V = 1.5 x 10" 5 where Kf = 4.8 x 10" 9 (left), and 
Wp F = -0.44, V = 0.075 where Kf = 7.6 x lO" 10 (right); 
together with the exponential fit described in the text. 
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where 5j z — nJ z p/2 is the phase shift resulting from </ z , 
and we have employed the bosonization cutoff scheme 
here. For small values of both J z and Jj_, the above 
equations can be expanded to yield: 

c = 4( J^pf - A{J zP f . (37) 

For J± — > the value of c thus depends only on J z as 
anticipated; in this regime c < 0. The advantage of using 
a parametrization of the single-impurity RG flow via c is 
that it allows to cover the trajectories with \J±\ > \Jg\ 
as well, i.e., the trajectories above the isotropic line in 
Fig. here c> 0. 

With the parameter mapping described at the begin- 
ning of this section we have all parameters at hand and 
can determine the value of c from V and W. This al- 
lows to re-plot the phase diagram in a plane spanned by 
/K z r and c - this is shown in the inset of Fig. [JJ In 
particular, we can now add data points for J z < to 
the phase boundary plot, as those are characterized by a 
finite J± in the limit -C D. 

As mentioned above, some NRG results show a rela- 
tively strong dependence on the NRG discretization pa- 
rameter A. Fig. [7| shows the phase diagram for A = 2; 
results for other A values are similar, but the T^/K? 
values can differ by 50% or more. Therefore, we have per- 
formed an extrapolation to A — ► 1 for a few important 
quantities. A sample extrapolation is shown in Fig. |S| 
for the slope of the phase boundary near the Toulouse 
point, which was determined analytically in Sec. II I II - 
the extrapolated value of K Z W jV 1 ^ is consistent with 
the exact result in Eq. I|31|) . We have also looked at 
the maximum value of Tj}'/K z of the phase bound- 
ary occuring near J z = 0, this value extrapolates to 

max — 

0.11 ±0.03. 



FIG. 7: Phase diagram of the generalized single-impurity An- 
derson model 11 1ft deduced from NRG calculations for NRG 
discretization parameter A = 2. The vertical dashed line 
shows the Toulouse point of the individual Kondo impurities. 
Small values of V have been used to reach the universal regime 
<C D. Precise values of T% have been determined via 
the specific heat coefficient 7, see text. The upper horizontal 
axis shows the corresponding values of J z in the bosonization 
cutoff scheme. The error bar shows the typical uncertainty 
in the numerical determination of /Kf arising from the 
fits of both 7 and Kf . The inset shows the same data for 
Tk /Kf, now plotted as function of the RG invariant c of 
the single-impurity model (I36H - this plot covers the range of 
positive as well as negative J z (here c > 0). The lines are 
guide to the eye only. 
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FIG. 8: A dependence of the slope of the phase boundary 
near the Toulouse point. (A is the NRG parameter defining 
the logarithmic discretization of the conduction band.) The 
dashed line is a linear fit. Each data point involves an ex- 
trapolation of the numerical results at finite negative W to 
W — > 0. A rather strong A dependence can be observed, 
however, the extrapolated value appears consistent with the 
analytical result 11311 . 
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V. FLOW EQUATIONS 

In this section, we consider a different approach 3 ^ to 
our original model of Ising-coupled Kondo impurities, 
which is based on the method of flow equations. 35 The 
general idea is an approximate diagonalization of each of 
the two Kondo impurities - the result is a resonant-level- 
type effective model which captures the Kondo physics in 
terms of a non-trivial renormalized hybridization. Tak- 
ing the two Ising-coupled impurities together, we will 
again arrive at an effective Anderson model. Away from 
the Toulouse line one finds an additional weak density- 
density interaction. However, in contrast to Sec. IIII Dl 
where we were forced to treat a similar interaction non- 
perturbatively to take into account x-ray edge singulari- 
ties (App.E}, this is not necessary here because this non- 
trivial physics is already contained in the fully renormal- 
ized couplings which naturally appear within the flow- 
equation approach. This approach allows for a system- 
atic expansion around the Toulouse line, and the effective 
Hamiltonian derived in this framework in fact describes 
the entire phase diagram of Fig. ^ consistently. 



A. Flow equation transformation 

The flow equation method was first applied to the 
Kondo model in Ref. l3fil where it was shown that it 
leads to an expansion around the Toulouse point. Its ba- 
sic idea is to perform a sequence of infinitesimal unitary 
transformations on a given many-particle Hamiltonian 
and thereby diagonalizing it. 35 The expansion parameter 
turns out to be Ao — 1 with Ao = \/2(l — J z p) = 1 — Wp 
(using the bosonization cutoff scheme); the Toulouse 
point corresponds to Ao = 1. 

Following Ref. l36l we construct unitary transforma- 
tions Ui t 2 (see App. [BJ) such that the single-impurity 
Hamiltonians H^ 2 from Eq. (|10|l become diagonal 

h[%~^ =U ia H^Ut a (38) 

up to higher-order terms in our expansion around the 
Toulouse line. For studying the Ising-coupled Kondo im- 
purities we therefore apply the combined unitary trans- 
formation U = Ui U 2 on @, H K = U H K U\ leading to 

H K = H[ K - dias) + H ( 2 K ~ dias) + K z SIS Z 2 (39) 

with SI 2 = U 12 Sl 2 U\ 2 . At the Toulouse point JSUJl 
is of course exactly equivalent to the Anderson impurity 
model (|13fl with W = and the same mapping as used in 
Sees. Ill (J and III Dl the unitary transformation £/i i2 just 
eliminates the hybridization coupling in the Anderson 
impurity model. 

It was shown in Ref. |33 that the flow equation ap- 
proach yields a resonant-level model (|13|l as an effective 



model for the Kondo impurity model also away from the 
Toulouse point 

if( Ri - eff ) =H [y j ] + Y,V k (d)* j {k)+Yi.c) , (40) 
k 

where the ^(k),^j(k) are the creation and annihila- 
tion operators for solitonic spin excitations in momen- 
tum space. However, this resonant-level model now 
has a nontrivial renormalized hybridization function, 

A(e) = Z k VZ5(e-e k ), with (i) A(0) = T«/^ 2 
and nearly constant in an energy interval of order 
around the Fermi energy (here ep = 0), and (ii) a non- 
trivial power-law behavior for larger energies. Further- 
more, it was shown in Ref. |33 that to leading order in 
an expansion around the Toulouse line one can identify 
Sj = d ] jdj — 1/2. The effective model for our system of 
Ising-coupled Kondo impurities is therefore an Anderson 
impurity model with a hybridization function of order 
the single-impurity Kondo scale: 

H (A-eS) = H {* a ]+Y,V k (49 t ,{k)+h.c.)+K,n dr n dl . 

k.a 

( 41 ) 

The main feature of the flow equation method is therefore 
to eliminate the large coupling W in (|14|l by renormaliz- 
ing the hybridization of the Anderson model. 

However, since the flow equation transformation is an 
expansion in the distance (Ao — 1) to the Toulouse line, we 
need to be careful in the transformation of S z : The trans- 
formed S z is multiplied by a possibly large parameter K z , 
so that an error of order (Ao — 1) in the expansion becomes 
multiplied by K z , leading to additional interaction terms 
in H41fl that can be larger than the hybridization energy 
scale T^p iSL It are precisely these additional interactions 
that drive the Kosterlitz-Thouless transition between the 
Fermi-liquid phase and the frozen mini-domain phase for 

\k z (\ -i)\>t£\ 

B. Corrections to the transformation of S z 

In App. [B] the flow equation solution for the single 
impurity Kondo model H^ 2 in a magnetic field h is dis- 
cussed with a careful analysis of terms of order (A — 1). 
Here h is the effective exchange field due to the second 
spin, to be described below. The transformed operator 
takes the following form 

Sj=\J dxdx' d(x)d*(x')[^ j (x),^ j {x')] 

+\{\ Q -l)f{h)d x fo{0) (42) 

plus irrelevant terms (containing e.g. higher derivatives 
of the bosonic field) and plus higher order terms of or- 
der (Ao — l) 2 . The first term on the right-hand side of 
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Eq. (|42|l can be interpreted as the result of integrating 
out the hybridization term in (|40|) , while the second term 
is a correction term not contained in the original solution 
in Ref. |3(| 4>j (0) denotes the bosonic spin-density field 
<j)j (x) without the Fourier components for energies larger 

than 0(Tft) (with respect to low-energy properties one 
does not need to distinguish these fields). Properties of 
the dimensionless function f(h) are derived in App. [Bj 



in particular f(h) = v F /\h\ + 0(h,- 2 ) for \h\ > . 

In the coupled system (|41[l we can approximate the ef- 
fect of one spin on the other as a static magnetic field of 
strength h = ±K z /2 close to the transition. This approx- 
imation becomes asymptotically exact as one approaches 
the transition since the spin dynamics becomes slower 
and slower. 



I 

We arrive at the following Hamiltonian describing the coupled Kondo impurities in the vicinity of the transition 
line. 

H (A- eS ) a H^^ + J^Vk (49 a (k)+h.c.)+K s n dt fi dl + (A - l)^f(K z /2) (d x ^(0)n dl + n dT x ^(O)) (43) 

k,(T 

up to corrections of order (Aq — l) 2 . 



The Kosterlitz— Thouless transition 



Let us now focus on the case K z ^> Tjp that is relevant for studying the phase transition in the vicinity of the 
Toulouse line. Using f(h) « vf/\Ii\ we rewrite the Hamiltonian l|43|) as 

H (A-eS) = H<j ^ + Y^V k (4%{k) + h.c.) + K z n d ^n dl 



-(Ao - l)v F ({d x fa(Q) + dMQ))~(n dl + n dl ) - (0«0 T (O) - d x ^(0))^(n d ^ - n dl )) . (44) 



r 



Similar to the analysis in Sec. IIII CI the term propor- 
tional to n d ^ + n d i is frozen out and can be ignored, 
while the term proportional to n d ^ — n d ± leads to a spin- 
spin interaction. Since K z is the largest energy scale in 
(|44|l with its renormalized parameters, we can map the 
Hamiltonian onto an anisotropic Kondo model using a 
Schrieffer- Wolff transformation (like in Sec. 1111 (J|) and 
again arrive at Ij23(l 



J2 jfs 11 



'^(0X^(0) , (45) 



na/3 



where now the Kondo couplings for scattering processes 
in the vicinity of the Fermi surface contain the renormal- 
ized parameters p F V 2 — /wtt 2 



PFJf 



4p F V 2 
K, 



p F J\ 



off 



Ap F V 2 
~K, 



PF 



J (nl) (46) 



with p F J^ — Ao — 1. We stress that here it was not 
necessary to use the generalized Schrieffer- Wolff trans- 
formation derived in App. El as the parameters in 144(1 
are already renormalized due to the flow equation proce- 
dure and the interactions oc Ao — 1 are only effective at 
low energies. The additional spin-spin interaction J*- nl ) 
is ferromagnetic for couplings to the right-hand side of 
the Toulouse line Ao > 1. This leads to a critical cou- 
pling for the Kosterlitz-Thouless transition to the frozen 



mini-domain phase 



k: 



8p F V 2 
A -l 



8V 2 
-W 



or using 



Eq. (EH 



PfK 



8 T. 



(i) 



A' 



-W 



= 1.964- 



T. 



(i) 



A" 



-w 



(47) 



(48) 



in exact agreement with the NRG results (Fig. |HJ and 
the strong-coupling analysis Eq. in Sec. IIII Dl 

Since the flow equation approach leads to a renormal- 
ized effective Hamiltonian, one can also use it to derive 
the entire phase diagram like in Fig. [7] If one neglects 
the same higher order terms (A — l) 2 as before, one finds 
the following result for the critical coupling 



PfK 



8 T, 



(i) 



K 



-w 



(1-A'(0)) 



(49) 



which simply results from a Schrieffer- Wolff transforma- 
tion of an Anderson model with an on-site repulsion K z 
and a hybridization A(e) which enters as . Here 



A(w) = ~P f de 



A(e) 



(50) 
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FIG. 9: Phase diagram of the generalized single-impurity An- 
derson model 11411 . deduced using the flow equation method. 
Notation is as in Fig. 



follows from the effective hybridization function of the 
flow equation approach—. The factor (1 — A'(0)) enters 
in 149|) because it generalizes the relation between the 
Kondo temperature defined in (|32|l and the renormalized 
hybridization function within the effective resonant level 
model 37 



ttA(O) = — (l-A'(Q)) 
w 



(51) 



The results are depicted in Fig. The maximum value 
of the phase boundary occuring near J z — is given by 
(T^ /K z )max = 0.126, which agrees with the extrapo- 
lated NRG value 0.11 ± 0.03 from Sect. EH 



VI. SYMMETRIES AND PERTURBATIONS 

To what extent do the results presented in the previous 
sections depend on the details of the models under con- 
sideration? To answer this question we will investigate 
whether and how (small) perturbations of © will qual- 
itatively change the physics. Fermi-liquid phases with 
vanishing residual entropy 5*0 are stable against small 
perturbations. This is not necessarily the case for our 
"frozen mini-domain" characterized by So = In 2. The 
existence of this So = In 2 phase is a fundamental fea- 
ture of our model J3J)- The necessary conditions for its 
stability will be discussed in what follows. 

Firstly, let us consider the effect of a magnetic field 
in z-direction acting on the impurity spins. A staggered 
magnetic field, h s (S^ — Sf ), will directly destroy the de- 
generacy of the two configurations | fi), | If). How- 
ever, in the limit K z -^ooa homogeneous magnetic field 
h(Sf + Sf) will not destroy the Sq = In 2 phase. It is 
interesting how these terms modify the generalized An- 
derson model (|14fl . The magnetic field h results in a term 
h d\ r d cr which breaks particle-hole symmetry in the 



generalized Anderson model. It therefore modifies only 
the position of the phase boundary. The staggered mag- 
netic field h s , however, will lead to a term h s J2 a CTC 4-^o- 
which corresponds to a (pseudo-) magnetic field acting on 
the pseudospin of the Anderson model. Only the stag- 
gered magnetic field is a relevant perturbation which de- 
stroys the In 2 phase. 

Apart from these magnetic fields in z-direction there 
are other relevant terms which lift the twofold degener- 
acy, which are of the forms: 



1,2, 



$1 S 2 , 



(52) 
(53) 
(54) 
(55) 



and their hermitian conjugates. It turns out that all these 
operators are forbidden if we impose the following two 
symmetry conditions: The model should be invariant un- 
der the two separate spin rotations of each impurity and 
its electronic bath about an angle of n, i.e., under the 
transformation 



U, 



(56) 



with j = 1, 2. Ij is the z-component of spin of system j, 
If 



Sf + E 



tt rotation symmetries, U 



k C kaj2 a af3 C kl3j 



J ' 



. In the presence of these 

the terms JSJ, 

and (|55|l are absent and the frozen mini-domain phase 
survives. The quantum phase transition from the frozen 
mini-domain with residual entropy In 2 to the phase of 
Kondo screened impurities therefore just relies (in the 
absence of a staggered magnetic field) on the symmetries 
U\ and U 2 - 

The model (01 considered in this paper possesses by 
construction symmetries beyond Uj . They are not neces- 
sary for the stability of the So = In 2 phase. For example, 
the two baths are assumed to have the same Kondo cou- 
pling J n . This parity symmetry can be relaxed without 
destroying the frozen mini-domain phase. Furthermore, 
the z-component of spin of each system, J? , is conserved 
in our model since we chose J x — J y — Jj_. This sym- 
metry can also be perturbed without lifting the twofold 
degeneracy. Moreover, the frozen mini-domain phase is 
stable against breaking of particle-hole symmetry which 
we implicitly assumed in the bosonization treatment by 
linearizing the dispersion relation of the conduction elec- 
trons. In all these situations, we therefore expect that all 
of the qualitative results, i.e., the structure of the phase 
diagram and the nature of the quantum phase transition, 
are not affected. 

However, any perturbation which breaks either U\ or 
U2 (or both) will generically generate one of the relevant 
couplings i|52H55|) which all destroy the In 2 phase. In 
the following we briefly discuss two such cases which are 
likely to occur in experimental realizations (a third case, 
corresponding to l|54|) is studied in Sec. IVII A"|l . 
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First, consider a situation where a small spin- flip cou- 
pling H53[) is added on top of the large Ising interaction 
of the spins, 

&H& = K ± (S?S% + S\S V 2 ) . (57) 

In realizations of our model based on spins and strongly 
anisotropic spin-orbit interactions, such a term will al- 
ways be present. A small K± will immediately lead to 
a tunneling between the two states of our mini-domain: 
their degeneracy is lifted, the two spins form a singlet and 
the In 2 residual entropy is quenched completely. Two- 
impurity Kondo models with K± = K z have b een wid ely 
studied. 12 ! 13 ' 14 ! 15 ' 16 !^! 18 ! 19 As argued in Refs. Il6ll7l the 
resulting phase diagram depends on the presence or ab- 
sence of particle-hole symmetry (which does, however, 
not modify the phase diagram for K± = as pointed out 
above). In the absence of particle-hole symmetry, the 
phase transition is replaced by a smooth crossover. How- 
ever, in the presence of particle-hole symmetry, the scat- 
tering phase shifts of the electrons can only take the val- 
ues or 7r/2. As the Kondo-screened phase and the inter- 
impurity singlet phase have different phase shifts, there 
has to be a phase transition in between. This transition 
is not of Kosterlitz-Thouless type, but characterizediS^ 
by a residual entropy In y/2. Nevertheless, this transition 
will merge with ours in the limit K± — > 0, as an infinites- 
imal K± does not affect the Kondo-screened phase but 
leads immediately to the formation of an inter-impurity 
singlet in the frozen mini-domain phase. 

A second interesting case is a situation where the two 
Fermi seas are coupled, e.g., by a tunneling between the 
two leads 

8B*r^ = E : c ! Ql ^ Q2 : + h-c.) . (58) 

k,k',a 

While this term is not relevant by power counting, it 
will induce an RKKY interaction between the spins and 
therefore generate the relevant coupling l(5S|) or (|5?|) . 
As such a term also breaks particle-hole symmetry, the 
quantum phase transition will be replaced by a smooth 
crossover. 



VII. TRANSPORT 

In this section we illustrate how the phase diagram 
and, more importantly, the corresponding quantum 
phase transition can be revealed in transport experi- 
ments. We shall discuss two experimental setups: (A) 
Capacitively coupled quantum dots, where the charge de- 
grees of freedom play the role of (pseudo)spinSfS* 9 " are a 
promising realization of our model. By adding a small 
inter-dot tunneling term, we obtain a characteristic zero- 
bias anomaly. (B) If the Ising coupling is realized be- 
tween real spin degrees of freedom, then we shall show 
that a transport experiment can reveal a universal frac- 
tional critical conductance at the phase transition which 



is related to the universal jump of the superfluid den- 
sity at the Kosterlitz-Thouless transition of superfluid 
thin films. Using quantum dots this situation is difficult 
to achieve, as a transverse spin coupling will always be 
present in the experiment. Nevertheless, the following 
proposals highlight the non-trivial effects of a Kosterlitz- 
Thouless transition with a bath of solitonic particles onto 
the original electrons. 

A. Zero-bias anomaly of capacitively coupled 
quantum dots 

In realization of our model JHJl using charge states of 
capacitively coupled quantum dots^ the conductance dis- 
cussed above cannot be easily measured. We therefore 
propose another experiment sketched in Fig. llUI We con- 
sider two large quantum dots, each coupled to a (single- 
channel) lead. The Coulomb interaction and an appropri- 
ately choosen gate voltage ensure that two charge states 
on each dot are degenerate, and that all other charge 
states have higher energies. These two charge states in 
each dot take over the role of the two spins as explained 
in more detail by Matveev>& The spin up and down states 
of the conduction electrons in our model correspond to 
electrons sitting either in the leads or on the dot, where 
we assume that the level spacing on this dot is small com- 
pared to temperature. Using this mapping, a capacitive 
coupling of the two dots directly corresponds to an Ising 
coupling The physical spin in such a system would 
translate to an extra channel index in our model. For 
simplicity we will, however, consider a situation where 
either strong spin-orbit scattering mixes those channels 
or where the spin is quenched by a strong magnetic field 
- in both cases we effectively deal with spinless fermions 
and thus with a single-channel model. 

We now consider a situation where the two dots are 
coupled by weak tunneling t in addition to the large inter- 
dot capacity. This tunneling takes the form 

H tun = <5+5 2 -*{ 1 (0)* i2 (0)+h.c. (59) 

We assume that the tunneling is sufficiently weak that 
it can be treated perturbatively in the experimentally 
relevant temperature range. This is precisely the situa- 
tion which was also considered by Andrei et alM Note 
that the approximation to consider only tunneling into 
^ai {% = 0) in Eq. (|59|l is only valid if the tunnel contact 
between the dots is sufficiently close to the lead contact 
(see Fig.UHl. 

We calculate the conductance in perturbation theory 
in the interdot tunneling t starting from the Kubo for- 
mula. The current through the link between the dots is 
then given by j = iS' 1 + 5 2 "(i*j 1 (0)* i2 (0)) + h.c. and the 
T dependence of the current-current correlator can be 
obtained from simple power counting. 

We first consider the "frozen mini-domain" phase. Fol- 
lowing the arguments given in Sec. IIHI the dimension 
of the tunneling term (or equivalently of the current 
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operator) with respect to this fixed point is given by 
dim[ff tun ] = dimtf] = 1 - (f ) 2 - (1 - 2 4f. Therefore the 
current-current correlator decays in time as t _2 ( diro b'] _1 ) ) 
and we obtain for the conductance 

G(T) ~ t 2 T- 2dim ^ = <2 r -4f(i-M) _ (6Q) 

This divergence of the conductance arises because the 
tunneling is a relevant perturbation which will finally 
destroy the "frozen mini-domain" phase and quench its 
residual entropy In 2 below some small energy scale. 
Eq. H60fl is therefore only valid for sufficiently small t 
when this scale is smaller than T. Furthermore, a finite 
domain-flip rate induced by (|18fl is required to obtain a 
finite current. Above we implicitly assumed that t is so 
small that the size of the current is solely determined by 
the smallest bottleneck for charge transport given by t. 

At finite voltage V ^> T, T in l|6UI) can be replaced by 
V and we expect a zero-bias anomaly characterized by a 
pronounced peak in the conductance: 

G(V)~ |y|- 4 ¥(i-¥). (61) 

Upon approaching the quantum phase transition, the di- 
vergence increases and at the KT transition takes the 
universal form 

G cr (T) ~ T- 2 ^" 1 ) «T-°- 83 , (62) 
G cr {V) ~ ~ |y|-o.83 (63) 

up to logarithmic corrections. 

In the Kondo-screened phase, we can calculate the 
qualitative behavior of the current-current correlator at 
the point in the phase diagram where K z — and 
the dots decouple. The current-current correlator then 
can be decomposed into two correlators of the form 
5+ (i)*j.(f),Sr (0)^(0)^ which decay asymptotically 

as 1/t. This can be seen if one identifies this correlator 
with the conduction electron T matrix (see Ref. jijj and 
references therein) which is characterized by a constant 
spectral density for low energies. The conductance there- 
fore approaches a constant for temperatures and voltages 
well below the characteristic crossover temperature T* to 
the Kondo-screened phase: 

G(V) ps G(T) ps const. (64) 

In Fig. I1UI we show schematically the nonlinear conduc- 
tance as a function of V in the vicinity of the quantum 
phase transition. 

In contrast to Eq. (JHOJ and Eq. ffify. N. Andrei et 
alMi obtained an exponentially small conductance in the 
"frozen mini-domain" phase and G ~ T A in the singlet 
phase, with which we disagree. 

B. Universal conductance of Ising-coupled 
quantum dots 

What is the most characteristic signature of the 
Kosterlitz-Thouless quantum phase transition which we 




FIG. 10: Left panel: Experimental setup to measure the tun- 
neling conductance between two capacitively coupled quan- 
tum dots. Right panel: Schematic plot of the zero bias 
anomaly of the conductance at T = 0. In the "frozen mini- 
domain" phase, S < 5c, the conductance diverges algebraically 
according to Eqn. 1611 . At the quantum phase transition, 
5 = 8 C , the exponent takes the universal value — 2(v2 — 1) 
according to Eq. Il(i2i . In the Kondo screened phase, 5 > 8 C , 
the conductance is finite for V — > 0. 



found in the previous sections? The most famous exam- 
ple of a Kosterlitz-Thouless transition is probably the 
vortex binding-unbinding transition in supcrfluid 4 He 
films. From the Kosterlitz-Thouless theory follows the 
prediction of an universal jump in the superfluid density 
upon passing through the transition^ 

Interestingly, the analogue of the superfluid density in 
our model is the scattering phase shift 6 of the conduction 
electrons, and the arguments for an universal jump in the 
superfluid density carry over to an universal jump in i5. 
This can be seen by considering the RG flow diagram 
Fig. |21 In the "frozen mini-domain" phase the system 
flows towards a line of fixed points which is naturally 
characterized by the dimension of the leading irrelevant 
operator, i.e. the domain flip l|18fl . or, equivalently, ac- 
cording to H21|) by the phase shift S of the conduction 
electrons. Upon approaching the quantum phase tran- 
sition, the irrelevant domain flips become marginal and 

the phase shift increases and reaches 5t = § f ^ — ^TfJ 
at the phase boundary [see Eq. (|22|l ]. On the other side 
of the phase diagram, the system flows to the strong- 
coupling fixed point, where the Kondo spins are screened 
and the electrons acquire a phase shift ofw/2. Therefore 
the phase shift jumps across the transition from 5t to 
7r/2! This picture is expected to hold everywhere close 
to the phase boundary as long as no other phase transi- 
tion intervenes - that the latter does not happen is shown 
by our NRG calculations. 

This universal jump of the phase shift has direct exper- 
imental consequences. Consider the experimental setup 
sketched in Fig. ^] where the conductance through the 
left dot is measured. If Kondo screening prevails, the 
conductance for T — > will be given by the conductance 
quantum Go = 2e 2 / \2ttK). In the frozen mini-domain 
phase on the other side of the phase diagram, spin flips 
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FIG. 11: Left panel: Experimental setup to measure the 
conductance through a single quantum dot, which is Ising- 
coupled to a second dot. The couplings to the leads and be- 
tween the dots can be tuned using appropriate gate voltages. 
Right panel: At T = 0, the conductance (solid line) takes 
at the quantum phase transition the universal value G cr = 
Go cos 2 -^m, 1661 . Dashed line: schematic plot of the conduc- 
tance at finite T. Corrections to the T = result are loga- 
rithmic at the transition. The exponent 2d = — 2dim[H^ p ] is 
given by the dimension of the domain flip term (11811 . 



are completely supressed for T — > and therefore we can 
assume a static spin configuration to calculate G(T = 0). 
For such a potential scattering problem, the conductance 
is given by 



G(T = 0) = G sin 2 5 . 



(65) 



Directly at the quantum phase transition, the conduc- 
tance therefore takes the universal value 



G cr (T = 0) = Go sin 2 S T = G cos 2 



2V2 



0.2 Go . 



(66) 



and it jumps to the Kondo value Go upon entering the 
Kondo-screened phase. This universal fractional conduc- 
tance at our quantum phase transition is one of the re- 
markable results of this paper. 

It is interesting to compare this to the well-known re- 
sult for the usual Kondo effect, where the conductance 
jumps from to Go when the exchange coupling J is 
tuned from ferromagnetic to antiferromagnetic. Both in 
Sec. 1111 Dl and Sec. IV CI we mapped our model close to 
the quantum phase transition to such a Kondo model. 
The fermionic degrees of freedom in these Kondo mod- 
els [Eq. (|23|l or Eq. I|45l) ] are, however, complex solitonic 
excitations in terms of the original fermions. While the 
phase shifts of those solitons vanishes at the quantum 
phase transition, the phase shift of the physical electrons 
takes the fractional value St leading to a fractional con- 
ductance. Also in other systems which are described by a 
Kosterlitz-Thouless quantum transition in terms of soli- 
tons, a universal fractional conductance of similar origin 
can be expected at the transition. 

In Fig. ^] the zero-temperature conductance close to 
the phase transition is shown. At any finite tempera- 
tures, the jump in the conductance is strongly smeared as 



sketched schematically in the figure. The T-dependence 
at lowest temperature is determined by the dimension of 
the leading irrelevant operators. In the Kondo-screened 
phase leading corrections for T — > to the Kondo con- 
ductance G are of order (T/T*) 2 for T < T* , where T* 
is exponentially small close to the quantum phase transi- 
tion (see Fig.EJ). In the "frozen mini-domain" phase, cor- 
rections to ||BSJ| vanish as T" 2dim [ ff "ff P ] where dim[i^ p ] 
is defined in Eq. 1211 Directly at the quantum phase tran- 
sition the exponent vanishes and leading corrections to 
Eq. I|66l) are of the order 1/lnT and therefore rather 
large. 



VIII. SUMMARY 

We have investigated a model of two Ising-coupled 
Kondo impurities using strong-coupling expansion, nu- 
merical renormalization group calculations, and a trans- 
formation based on the method of flow equations. Those 
methods yield consistent results and allowed us to show 
the existence of a Kosterlitz-Thouless phase transition 
between a Fermi-liquid phase and a pseudospin dou- 
blet phase which corresponds to a "frozen mini-domain" . 
This transition can be tuned both by varying the Ising 
coupling between the impurities and by varying the 
anisotropy of the individual Kondo couplings. In par- 
ticular, at the Toulouse point of the individual Kondo 
impurities we could map the model exactly to an An- 
derson impurity model with a Fermi sea consisting of 
fermionic soliton excitations - in this situation no phase 
transition occurs, and the system is in the Fermi-liquid 
phase, where the impurity pseudospin is screened below 
a collective Kondo scale T*. For J z smaller than the 
Toulouse point value, large K z drives the system into 
the pseudospin doublet phase. 

The most promising way to realize our model is the sit- 
uation of capacitively coupled quantum dots where the 
impurity spins represent charge degrees of freedom on the 
dots. We have shown that a small additional tunneling 
between the dots gives rise to a zero-bias conductance 
anomaly with a universal fractional power-law occurring 
at the transition point. In addition, we have discussed a 
setup which is interesting on theoretical grounds, namely 
transport through one quantum dot of a pair of dots with 
a magnetic Ising coupling, where we have found a uni- 
versal fractional conductance through the device at the 
phase transition point. 

With an eye towards comparison with experiments we 
discuss the finite-temperature crossover behavior across 
the phase diagram (see also Fig. [SJ) If we fix the param- 
eters of the individual Kondo impurities, then varying 
K z corresponds to a vertical cut through the phase dia- 
gram in Fig. ^ the resulting finite-temperature behavior 
is sketched in Fig. IT21 

For small K z there is a single crossover at the single- 
impurity Kondo temperature Tg . This crossover splits 
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FIG. 12: Schematic phase diagram as a function of K z and T 
for fixed T^. 1 ' . For T — there is a quantum phase transition 
at K z = AJ r from a Fermi liquid with residual entropy So = 
to the "frozen mini-domain" phase with So = In 2. At T > 0, 
only smooth crossovers occur indicated by the dashed and 
dotted lines. At the dashed lines, the entropy S changes by 
In 2 (see also Fig. 0, while at the dotted line one obtains 
a crossover from a logarithmic to a power-law behavior in 
the leading corrections to S. Similar crossovers also occur in 
transport quantities. Below To a magnetic "mini-domain" is 
formed, while a Fermi liquid is recovered below T* which is 
exponentially small close to K^ r . 



into two when K z approaches values of order Tj^ 



then 

the described two-stage quenching of the entropy is ob- 
served. The upper crossover temperature, To, is associ- 
ated with the formation of the magnetic mini-domain, 
where relative fluctuations of the two impurity spins are 
frozen out. The lower crossover temperature is the col- 
lective Kondo scale T* below which the pseudospin of the 
mini-domain is screened. T* becomes exponentially small 
near K z r and vanishes for K z > K z r . For K z > K z r an- 
other crossover line appears which, however, has much 
weaker signatures, namely the character of the leading 
corrections to the entropy and other quantities changes, 
as is easily understood from the RG flow in Fig. [3J For 
large K z the entropy change from In 4 to In 2 occurs 
around T ~ K z , therefore T approaches K z in this limit. 

Interestingly, the different impurity degrees of free- 
dom can be re-interpreted: the flipping of the pseu- 
dospin while keeping the mini-domain intact apparently 
corresponds to pseudospin "phase" fluctuations, whereas 
breaking up the mini-domain is related to "amplitude" 
fluctuations of the pseudospin. Thus, in Fig. 1121 we 
encounter the situation that amplitude fluctuations are 
frozen out at a higher temperature T whereas phase fluc- 
tuations are quenched at the lower T*, in other words, the 
two impurity spins fluctuate independently for T > T 
whereas they fluctuate in a correlated fashion between 
T* < T < T . This physics is surprisingly similar to the 
behavior of lattice systems in low dimensions, with the 
difference that of course no true ordering can occur in 
the impurity model. 

In summary, the present two-impurity model shows 
remarkably rich behavior, which awaits realizations in 
mesoscopic devices. An interesting extension would be 
the two-channel case which is naturally met in capaci- 
tively coupled dots with spin-degenerate conduction elec- 
trons. 
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APPENDIX A: GENERALIZED 
SCHRIEFFER— WOLFF TRANSFORMATION 

In this appendix, we perform explicitely the mapping 
of the generalized Anderson model (|14f> to the Kondo 
Hamiltonian (|23fl for large K z . Due to the presence of the 
interaction W in i|14|) the usual Schrieffer- Wolff transfor- 
mation has to be generalized to take into account power- 
law renormalizations (of "x-ray edge" type) induced by 
W . We derive the mapping by investigating directly the 
properties of a perturbative expansion in the hybridiza- 
tion V for finite W within the Anderson model. Consider 
the generalized Anderson model in its bosonized version. 
We first eliminate the W term in i|14|) by the Emery- 
Kivelson transformation JjJJ with 7* = W p and obtain 



U^H A U^=^2H y> tr ]+K. 



H 



Lnt 



(Al) 



where W enters only the hybridization term 
V 



Hi nt = 



; \-^t e - l (i^(o) Fo+La j _ (A2) 



For large K z the d-level is only singly occupied. V in- 
duces virtual fluctuations to the the doubly occupied and 
empty state which are separated from the singly occupied 
states I t) , I I) by an energy K z /2. To derive the effective 
Kondo model consider the S-matrix with respect to this 
low-energy subspace: 



00 

Texp[-i J dtH int (t)] 



E 

n=0 



dt 2n ■ ■ ■ dtxiH iDt (t2n) ■ ■ ■ iH ia t(ti) (A3) 

t2»>-->tl 



Hi nt describes processes from the low-energy sector to 
high energies or back. Such virtual excitations are rare 
and exist only for a short time if K z is large. Therefore we 
can group them to pairs to obtain an effective interaction 
living in the low-energy Hilbert space, 
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*2m+l t 2m 

dt2m / dt2m-liHi nt (t2m)iHi nt (t2m-l) 



, r m+1 Hg*(T m ) . 

OO 



(A4) 



OO 

with H^ t {T m ) = -ifdt H int (T m + t/2)H int (T m - t/2) 



where we introduced the center-of-time and relative coordinates. Interactions between adjacent virtual excitations 
can be neglected to leading order for large K z . Introducing the spin notation, S z = \ ^ CT <jS a d a and S + = did,, to 
represent the two states of the low-energy Hilbert space, the above expression becomes 



V 2 f°° 

-i / dte 

2na J 



q i2S z o-(l-Wp)<p„(T m +t/2) e i2S z o-(l-Wp)<p< 7 (T m -t/2) 



a 

+ (S + F F^e'^ 1 '™ »<M T ™+*/ 2 ) i e *o-(l-w ' P )4>-<,(T m -t/2) 



-^int(-^m) 



The oscillating factor e lK » t / 2 guarantees that the virtual excitations are only short lived, and we can therefore 
expand the term in the bracket in the small time t. Introducing the spin field </> = ^= cr(f> a , applying the operator 

product expansion e iMa{t) e -i\<t>a{t') = (1 + i(t - t')/a)~ x2 + Xa(l + i(t - t') / a) 1- * 2 'd t '(j) a {t') + ... for the first term, 
integrating over t using J °° dte~ lKzt ^ 2 (l + it/a)~ a = ~i(aK z /2) a 2T(l — a)/ K z we obtain in leading order for large 



H '» = l^ , ™ (2 - (1 - HV,3, ("? i ) 



(i-ffp) -i 



dxm : +^L( S +e- i ^ 1 - w ^F}F, +h, 
K z 2ira \ 1 1 



r 



iA5) 



Before identifying the coupling constants of the effec- 
tive low-energy Hamiltonian two more steps are required. 
First, we have to re-adjust our UV cutoff from a to 
o-k ~ m the definition of our fields, as we effec- 

tively have integrated out short time differences of order 
l/K z . To this end we have to normal-order H?f t as only 
normal-ordered expressions are cutoff independent, 



iX<p 



27raV 
~L~J 



: e 



a K 



V L 



a 

a K 



(A6) 



where 4> denotes the fields defined with respect to the 
new cutoff aK- This effectively leads to the substitution 

(1-Wp) 2 -1 

(A7) 



AV 2 



W 2 



- 



K z 2na 



K z 2ira,K \ax , 

in the second term of (| A5|) . 

In a last step, we undo the Emery-Kivelson transfor- 
mation to obtain the Kondo Hamiltonian in its usual 
form with 



cir 



j: 



jf = 



w + 

AV 2 



4V 2 
a 

ax 



-cw 



aK- 



2 

{i-WpY 



(1-Wp) 2 



(A8) 
(A9) 



where c w = (1 - Wp)T(2 - (1 - Wp) 2 ). X-ray edge 
singularities induced by W have led to a power-law de- 
pendence of the effective couplings on K z . Note that the 
previous arguments fixed aK ~ l/Kz m <|A9|) only up to 
a prefactor of order 1 depending on W. However, this 
unknown prefactor approaches 1 close to the quantum 
phase transition where Wp — + 0. 



APPENDIX B: FLOW EQUATION 
TRANSFORMATION FOR THE SINGLE KONDO 
IMPURITY 



In this appendix we provide some details on the flow 
equation treatment of the single-impurity Kondo model, 
which was first presented in Ref. I'itii Here we will show 
how to extend the analysis of Ref. |3^ to take into ac- 
count the terms of order (A 2 — 1) that become important 
in our coupled system since they are multiplied by a pos- 
sibly large energy scale K z in 1)41(1. We refer the reader 
to Ref. 36 for the basic ideas of the approach and only 
present the main steps to keep our presentation here self- 
contained. 
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1. Transformation of the Hamiltonian 

The starting point for the flow equation approach is 
HI U| l with 7 chosen such that the longitudinal coupling is 
eliminated. This way we arrive at the initial Hamiltonian 
H(B = 0) for the flow equation approach 



H(B)=H [<f>] + J dxg(B-x) (V(X(B);x)S~ +h.c. 



with X(B = 0) = A = a/2 - J z /V2nv F and 



(Bl) 



/9 v X(B=0) 2 /2 j 

g(B = 0;x) = S(x)(^j (B2) 
Here V(X;x) are normal ordered vertex operators 



V(X;x) 



— i\<p(x) . 



(B3) 



During the course of the infinitesimal unitary transfor- 
mations 



dH(B) 
dB 



= [r,(B),H(B)} 



(B4) 



with the generator rj(B) from Ref. |3a 

t] = J &xri {1) {x) {V(X;x)S- - h.c.) (B5) 
+ / dxdx'r)W{x,x') [V{X;x),V(-X;x')] 



the interaction g(B;x) in l)Bljl becomes more and more 
nonlocal. With each infinitesimal step of the transforma- 
tion one also generates a new interaction term in i|Bl|) 
with the structure 



vf S z / dx s(x) d x <f>(x) 



(B6) 



and a nonlocal function s(x) that depends on the cou- 
plings. The key step in Ref. |36| is that ljB6(l can again 
be eliminated by a unitary transformation of the Emery- 
Kivelson type 



(B7) 



U = : cxp iS z / dxs(x) 4>(x) 



We now analyze how the interaction term in 1|B1|) is transformed due to U, e.g. 



UV(X;y)S-tf =: cxp ( -- / dxs(x) <f>(x) 



e- lA0(s/) : : exp 



•- / dxs(x) cf)(x) 



(B8) 



In order to proceed we normal order all the exponentials, which can be done exactly since the commutator of the 
bosonic field is a c-number. This leads to 



UV(X ] y)S-U^:e X p^ T Y / ^ r ([^ k y + s(k)]b k ~[Xe^ +s(k)]b{)) : S~ 



(B9) 



with s(k) being the Fourier transform of s(x) from (|B6|I . The proportionality factor in l|B9|l leads to the non- 
perturbative renormalization of the coupling constant g(B; x) already obtained in Ref. 1361 Except for the local 
coupling at the beginning of the flow the exponential in <|B9|) cannot be exactly rewritten as a vertex operator. We 
use two approximations that give us the correct result up to quadratic terms in the deviation from the Toulouse 
line: i) We use the infrared limit s(0) instead of s(k) in (|B9|I . ii) We expand the exponential in a way that avoids 
IR-divergences and neglect higher order terms in the bosonic operators that lead to irrelevant couplings: 



: ex p (Vt S ^ ([Ae ~ lfe " + s(0)] K _ h - c - } ) 



: exp 



E 



-ka/2 



L t?0 ^ 



-([(A + s(0))e"*f 



= ■.V(X + s(Q);y)(l + J^J2 [ 

k>0 



-ka/2 

~7k" 



-e- ft »)a(0)]6 fc -h.c.)) : 
^) S (0) 6 fe -h.c.)+...) 



(BIO) 



Retaining only the first term on the right-hand side is eventually become fermionic. The second term can be 
the approximation used in Ref. l36t one obtains vertex understood as a correction term to this leading behavior 
operators with flowing scaling dimensions [A + s(0)] that due to the nonlocality of the interaction during the flow 
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equation procedure. It is this term that eventually leads 
to the correction term in (|42|l . 

The above procedure following from l|B8(l has to be 
repeated iteratively throughout the flow, leading to 

x((l-e-^)6 fe -h.c.)) :S~ 
= :&(y)(l + (l-\ )(-im+i$(y))) : S~ 
= : ¥{y) (l + (1 - A ) iyd v m) : S~ (Bll) 

plus irrelevant terms with higher order derivatives of the 
bosonic field. Here 4>{y) denotes the bosonic spin-density 
field 4>(y) without the Fourier components for energies 
larger than O(Tr-) since the term proportional to (1 — 
Ao) is generated successively during the flow equation 
procedure. Putting everything together, the Hamiltonian 
H{B) from (|fjl() acquires a new term of order (1 — Ao) 
during the flow that has been neglected in Ref . |3|| It can 
be viewed as an assisted hopping term that is marginal 
as opposed to the leading order hopping term that is a 
relevant operator. The new term can be eliminated by 
including an additional term with the structure 

J dx rf® (as) ( : V(X; x)d x 4>{0) : S~ - h.c.) (B12) 

and a suitable coefficient function rj^ (x) into the gener- 
ator (|B6|> . One can verify that this does not modify the 



previous flow equations for the Hamiltonian in linear or- 
der in ?y( 3 ) (x) (essentially since the assisted hopping term 
is marginal as opposed to the relevant hopping term that 
generates the flow equations in leading order). There- 
fore we can neglect these extra terms in the flow of the 
Hamiltonian when we want to retain terms up to linear 
order in (Ao — 1). 



2. Transformation of the impurity spin operator 

However, for the transformation of S z one needs to 
be more careful since S z can be multiplied by a large 
exchange field h — ±K z /2 due to the coupling to the 
second spin. This can be much larger than the Kondo 
scale close to the transition. In order to study the trans- 
formation of S z we follow the same route as in Ref. |3^ 
by using the identity S z = {S + ,S~]/2 and evaluating 
the transformed S + (S~ then follows as its hermitean 
conjugate). One needs to study the additional effect of 
(|B12|) on S + and finds the following expression in the 
low-energy limit: 

S + = <J Z J dyd(y) : ¥(y) (l + (1 - A )i9^(0)) : 

(B13) 

with [to linear order in (Ao — 1)] the same coefficients d(y) 
as in Ref. |3(| This leads to 



S z = [S + ,§-]/2 = - J dxdx' d(x)d*(x')[^(x),^{x')} + ^(l- \ ) J dxdx' d(x)d*(x') 

x ( - ix'[tft(x), : *(x')d x 4(0) :] + ix[: *Hx)d x 0(O) :, *(x')]) + O((A - l) 2 ) . (B14) 

Since we are interested in an analysis in the vicinity of the Toulouse line we only keep terms up to linear order in 
(Ao — 1). The term of order (Ao — 1) consists of two fermionic operators and a spatial derivative of the bosonic field. 
If we subtract the contractions with respect to the ground state the remaining normal ordered operator will therefore 
lead to an irrelevant coupling in the coupled Hamiltonians (|39|l . However, we need to retain the contractions: 

S z = i J dxdx' d(x)d*(x')[¥(x),y(x')] + ^(l- \ o )id x 0(O) J dxdx 1 d{x)d*{x') (x-x')([^(x),y(x')}) 
+O((A - l) 2 ) + irrelevant 

= 1 [dxdx 1 d(x)d*{x')[^(x),^(x')} + ^(Xo-l)f(h) d x 0(O) + O((A - l) 2 ) + irrelevant (B15) 
Z J 2 

with 

f(h) = J dkdk' [d k d k dy+d k d k ,d k ,) . (B16) 

Here d k denotes the Fourier transform of d(x). The expectation value ([^L ^fe']) has to be evaluated in the ground 
state of the Hamiltonian H ( pot ) that is obtained from the resonant level model Hamiltonian plus magnetic field h S z 
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after the above unitary transformation. Since S z "decays" into fermion operators under this transformation according 
to HB14|) . this Hamiltonian is given as 



/?&>«*) = £ e k *l* k + £ h d k dy , 



(B17) 



k.k' 



i.e. this is just a potential scattering model with a separable potential V kk i — hdkd k >. The retarded Green's function 
can be calculated in closed form 



leading to 



G kk ,{e + ) = 



5kk' 



h d k d k ' 



£fe (e+ - e fe )(e+ - e fe /) 



f(h) = J dkdk 1 (d k d k d k , +d k d k ,d k ,)([vl* k ,]) 

= -~Im.J dkdk' (d k d k d k , + d k d k >d k >) yj &eG kk >(e + )-J deG kk >{e + ) 



-Im 



-h — Im 

7T 



de 



oc 




de dk 



d k d\ 



de 



de dk 



dkd 2 k 
e+ - e k 

I 



dk 



, 4 



^9 £+-£, 



(B18) 



(B19) 



One easily shows that the impurity orbital Green's func- where p^ d °' ) (e) is the impurity orbital density of states 
T dd 



tion G, e t in the resonant level model 



One notices that the impurity orbital Green's function 
in the resonant level model with nonvanishing impurity 
H (RLM) _ ek^l^k + ^2 V( d ^k + ^k d ) ( B2 °) orbital energy e d d^d can be written as 



is given by 



G 



dd 



Using 



,+ , - VF d ^G dd (e j 



(B21) 



(B22) 



one can reexpress (|B19(1 as 
f(h) = 2 WF p^= 0) (0) 



x4r°V) 



which leads to 



i-e d 4r 0) (^) 



(B24) 



/W = 2«F^ £d - 0) (0) 



-Uf h — Im 

7T 



de 



de 



x(deG%r 0) (z + ))G d T h \t + )- (B25) 
This expression can be easily worked out in various limits 



2wv f /T ( k ] for/i = 
u F /|/i| for > T^ 1 



(i) 



(B26) 



i-/*4r°V) 



and a smooth crossover in between (here w = 0.4128 is 
(B23) the Wilson number). 
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